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Abstract 

Distance geometry is the study of the arrangements of points in space using only the mutual dis- 
tances between them. The basic idea in this letter is to use distance geometry for thermodynamics 
studies of small clusters in the microcanonical ensemble. There are constraints on these distances, 
which are shown to explain some characteristic features of the caloric curve in very small clusters 
containing 3 or 4 atoms. We anticipate that this approach could give a novel insight into the phase 
transitions in larger clusters as well. During these studies, we have established a very general 
and rather simple result for the Jacobian determinant of the change of variables from Cartesian 
coordinates to mutual distances, which is of wide applicability in the A^-body problem. 

PACS numbers: 36.40.Ei, 64.70.dm, 05.20.Gg, 82.60.Qr 
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It is well known that most of the tractable problems of physics are such because there 
exist a "good" choice of variables, which drastically simplifies equations to be solved. Such a 
good change of variable has not been found yet for the A^-body problem. Actually, for com- 
putational purposes, it has been agreed that the Cartesian coordinates are the less expensive 
coordinates because of the simplicity of the equations of motions when written in terms of 
them, at least in classical mechanics. On the other hand, those coordinates introduce unnec- 
essary invariance (the euclidean group), while physics describes mainly interactions between 
systems characterized by their relative position. In this respect, mutual distances should be 
enough to describe all physical processes. Another reason for considering mutual distances 
as variables is the simplicity of the potential energy when written in terms of them. It might 
thus be worth the try to go as far as possible with these variables in the studies of model 
systems described by pair potential. The work in this letter is a start in this direction and 
primarily of exploratory nature. It aims at bringing attention to some remarkable results 
(the simplicity of the Jacobian of the transformation from Cartesian coordinates to mutual 
distances) and to the ease of representation given by these variables. 

Mutual distances have been considered by Lagrange in his study of the 3-body problem. 
Several formulas for expressing various quantities of the system as a function of distances 
alone are known. The basic reference for these is Blumenthal's book A more modern 
account is given in chapter 4 of . Note that distance geometry is an invaluable tool in the 
study of protein conformation, where the problem is to determine all the distances between 
a set of atoms, knowing just a few of them, as given by NMR or X-ray studies 3, J]. 

Several new directions can be explored and some results will be presented in forthcoming 
publications. In this letter, this change of variable is applied to the computation of the 
configurational density of states (CDOS) of a small cluster of A^ atoms interacting by a pair 
potential, in a D dimensional space for the special cases N = DotN = D + 1. In both 
cases, the number of distances is equal to the number s of independent variables specifying a 
configuration: s = DN — D{D + l)/2. Note that for A^ > -D + 1, one has to choose a subset of 
the distances as internal variables, which introduces technical complexities not usefuU for the 
purpose of this letter. The general case will be presented in a forthcoming publication [5j . The 
CDOS is the fundamental quantity in classical thermodynamics. Its computation involves 
an s-dimensional integral. We show below that, when mutual distances are used as variables 
for an A^ = D-atom clusters, this s-dimensional integral can be approximately reduced to s 
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convolution products. This reduction is even exact at low energy. In any case, it is much 
easier to compute than an s-dimensional integral. 

The N = D = 3 case is especially interesting because it is possible to draw a represen- 
tation of the constant energy surface in the distance space. The accidents of the density 
of states as a function of energy (or rather of its logarithm, the entropy S{E)) can be vi- 



sualized and given a simple interpretation on this drawing. It is we 



1 known that accidents 



of S{E) are related to so-called phase transitions in clusters [g], Q, s], I^. The solid-liquid 
transition (actually a phenomenon akin to this transition) has been observed both in simu- 
H] and in experiments l3, Q, Q, 15, 16, 17|. We anticipate that the simple 
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view given in this letter could be used in understanding the solid-liquid transition in larger 
clusters as well. 

This letter is organized as follows. First, we recall the formulas for the microcanonical 
thermodynamics of a set of atoms and the definition of the CDOS. Then, we make some 
considerations on the change of variables, and give the (very simple in our case) formula 
for the Jacobian determinant. We next give some properties of the space of distances. We 
then present drawings of the constant energy surface in this space in the case of = 3 and 
explain the qualitative features of the entropy curve. Finally, we show how the CDOS can 
be expressed as a convolution product and present a comparison with more conventional 
calculations. 

In the microcanonical ensemble, the system size (volume V, number of constituents A'^) 
and energy E are fixed. The fundamental function is the density of states Q{E,V,N). 
It exists for any finite system with a classical Hamiltonian if(P,X). Here, the notation 
X represents the set of Cartesian coordinates of the constituents and P the corresponding 
momenta. The Hamiltonian H is the sum of the kinetic energy T(P) and the potential 
energy Ep{X), so that the density of states is actually the convolution of the kinetic density 
of states and the configurational density of states Qc{E, V, N) 

n^{E,V,N) = ^[ d'^X. (1) 

Since the kinetic density of state can be computed analytically, the configurational density of 
state is the quantity we focus on. The configurational entropy Sc{E) is defined by Boltzmann 
formula 

Sc{E, V, N) = lnfie(^, V, iV), (2) 



with kg the Boltzmann constant. In the remaining of the paper, the configurational entropy 
will be called just entropy for short. 

It can be shown [18] that, under mild conditions on the potential energy function Ep{X), 
the entropy is extensive in the thermodynamic limit, and that in this limit, the second 
derivative of the entropy is always negative. This is not true in a finite system: the entropy 
is not extensive because surface effects are important, and furthermore it can show portions 
of positive curvature, which are the accidents of the CDOS we mentioned above. 

We now describe shortly the change of variables from Cartesian coordinates to mutual 
distances. The Cartesian coordinates of an atom a, (a = 1 . . . A^), are noted x^, {i = 1 . . . D). 
Rather than the distances, it is more convenient to define the variables i/ab by 



Vab 



D 
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(3) 



When going from Cartesian coordinates to internal coordinates, one should define a "moving 
frame" (see for example Ref. 1^. Amazingly enough, the Jacobian determinant for the 
transformation does not depend on this choice. We suppose that the cluster is contained 
in a (hyper) spherical volume of radius R, with the center of mass fixed at the center of the 
sphere. This eliminates the need to integrate on the translation degrees of freedom. The 
integration on the rotation coordinates gives a constant irrelevant multiplicative factor, the 
volume of the group SO{D), noted Vrot- Eventually, it can be shown that the configurational 
volume element can be written 



/ JJ dxia = Hot J Yi 

^rotation •„ 



(4) 



l<a<6<Af 

with J' = 1 if N = D and J7 = 1 /V if N = D + 1, where V is the volume of the parallelotope 
built on the D + 1 points. This volume can be computed as the square root of a determinant 
involving the i/ab^ only |]|: 
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FIG. 1: (color on line) The distance space for = 3 atoms. The axes for the coordinates 
(yi2, yi3) ^23) ai'e shown by the intersection of two planes. An acceptable point can only be inside 
the cone. Furthermore, if the cluster is contained in a sphere, the accessible space is further limited 
by 3 planes, which interset on the (1,1,1) direction. On the figure, one of the planes has been 
made transparent in order to show the cone interior. 

Let us now think of the set of distances Uab as defining the Cartesian coordinates of a point 
in an abstract space which we call the 'Distance Space'. There are several conditions to 
be verified by this point. They amount to ensure that all the determinants similar to that 
of Eq. ([5]), for any subset of the N atoms, are positive. It can be shown (see e.g. Ref. [2I) 
that these conditions restricts the representative point in the distance space to be inside a 
convex cone. In other words, P is a space with a frontier. 

We have not yet included the condition that the atoms are confined inside a sphere of 
radius R. Actually, the distance from atom a to the center of mass can be expressed in 
terms of the y variables alone: 



Imposing that < amounts to adding the further constraints that the distance space 
is limited by a set of (hyper)planes. Those planes, together with the cone defined above 
define a finite volume in the distance space. Figure [T] shows a view of the cone and the 
planes in the case = 3. 

Now, What is the shape of the constant energy surface (CES) in distance space? We 
suppose that the potential is Ep = '^a^b'^iyab): where v{y) has the usual shape of a pair 
potential: one minimum at y = 1 with value -1, the limit when ?/ — ^ 00 is (that is the 
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dissociation energy is 1) and the limit when y ^ is +00. The CES in the case = 3 
is shown in the insets of Fig. ([2]). At low energy (inset (a), Ep = 0.2), the surface is 
close to a spheroid around the absolute minimum of the potential energy, which is yi2 = 
yi3 = 2/23 = 1/2. Notice that this surface is far from the limiting cone. At higher energy 
(inset (b), Ep = 0.5), the surface becomes elongated in the direction parallel to the (distance) 
axes. More precisely, it is easy to see that the intersection of the surface with the axis 
yi2 = yi3 = 1/2 for example, occurs for 7/23 such that f(?/23) = Ep. This value of j/23 
increases rapidly with increasing Ep. There is of course a symmetry between the three axes. 
On the other hand, at direction 45° from the axes, the maximum value of the y's is such that 
v{y) = Ep/2, which is much closer to the point (1/2, 1/2, 1/2). At an energy close to 1, the 
shape is that of 3 "tentacles" along the directions parallel to the axis, and the surface reaches 
the enclosing cone for Ep ~ 0.969 (inset (c)). This last value is specific of the Lennard- Jones 
potential; it is always very close to one for any potential. For E > 1 (inset (d)), the CES 
is open. However, the volume inside the cone is still finite. The configurational entropy is 
shown on the same figure. At low energy, Sc{E) oc ln{E), which is the usual behavior. At 
higher energy, the entropy departs from the logarithm, and would even diverge when E 1 
(the dashed curve). However, the limitation imposed by the cone results in a discontinuity 
of the slope when the CES meets the cone. It is important to understand that there is a 
balance between the tendency of the pair potential to allow a distance to become large, which 
results in a positive curvature of Sc{E), and the limitation of pure geometric nature due to 
the limiting cone, which bends in the negative direction. Such "bumps" of the entropy curve 
are characteristic of a phase transition in a finite system. We give here a new geometric 
view of what happens, at least in this simple system. 

Let us now compute the CDOS. The simplest case is N = D. When e < 1, the equation 
v{y) = e has two solutions, which we denote ?/max(e) and ?/min(e)- When e — > 1, ymax 00 
and disappears for e > 1, while ymin exists for any value of epsilon. To begin with, we forget 
about any constraint on the distances, and we suppose that N = D. It is easy to see that 
the CDOS can be written in this case 



Such an integral amounts to s — 1 convolutions. If the function uj is given on n points, 




(7) 



with Lo{e) = — [?/max(e) - ymin(e)]- 



(8) 
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FIG. 2: Configurational Entropy 5 as a function of energy E for a 3 atom cluster interacting via 
Lennard- Jones potential. The insets represent the CES in distance space (see text): (a) E = 0.2, 
(b) E = 0.5, (c) E = 0.969, (d) E = 1.2. The dashed curve is the entropy without constraint, the 
dotted curve is the result of a Monte-Carlo calculation, the plain curve is the approximation by a 
convolution product (see text). 

the integral can be computed in 0{n'^s) time, to be compared to the 0{n^) necessary for a 
general s-dimensional integral. Equation is valid as long as the energy is low enough for 
the CES to not intersect the limiting cone. Above that, an approximate density of states 
can be computed if one supersedes Eq. by 

a;(e) = ^[?/max(e) - l/min(e)] if l/max(e) < 1^, (9) 
= ^[ -Z/min(e)] if y„,ax(e) > 1^. (10) 

Y = {N — 1)/N — 2) is the value that zeroes the determinant ([5]) with all the y's equal to 
1/2 except one, which is equal to Y. The result of this approximation is shown on Fig. [2] for 

= 3 and Fig. [3] for = 8. Notice that the good agreement with the exact CDOS extends 
far beyond the expected range, which is E < 1. 

When N = D + 1, the integral in Eq. contains one more factor, the Jacobian determi- 
nant, so that it is no more a convolution. On the other hand, the CES develops essentially 
along directions parallel to the axes, hence all the y^s are close to 1/2. It can thus be included 
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FIG. 3: Configurational Entropy S" as a function of energy E for an 8 atom cluster interacting via 
Lennard-Jones potential. The dotted curve is the result of a Monte-Carlo calculation, the plain 
curve is the approximation by a convolution product. 
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FIG. 4: Configurational Entropy 5 as a function of energy E for a 3 atom cluster interacting via 
Lennard-Jones potential in dimension 2. The dotted curve is the result of a Monte-Carlo calculation, 
the plain curve is the approximation from Eq. (|lip . 



in the integral by changing one of the tu's with 

1 dy 

ma> 
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(11) 



V^bmax(e)] de V[ymm{e)] de 
where V{y) is the volume of Eq. (jS]) with all the yabS equal to 1/2, save for one, which is 
y. Of course, when y^ax > Y , only the term containing ?/min is retained as above. This 
approximation together with the exact Monte-Carlo result is shown on Fig. H] for 3 atoms 
in dimension 2. The agreement is still very good, although the approximation is somewhat 
crude. 

As a conclusion, mutual distances have been used as variables for the study of some A^- 
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body problems. It has been shown that it is possible to make efficient approximations for 
the calculation of classical densities of states and that considering the space of distances is 
a good help for viewing what happens at phase transitions in small systems. 

I thank Ms. Ch. Lauzeral for help on the initial stages of this project. Many calculations 
have been done or checked with the help of the Maple software from Maplesoft, a division 
of Waterloo Maple Inc. 
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